LONG TIME STABILITY OF A CLASSICAL EFFICIENT SCHEME 
FOR TWO DIMENSIONAL NAVIER STOKES EQUATIONS 
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Abstract. Wc prove that a popular classical implicit-explicit scheme for the 2D incompressible 
Navier-Stokes equations that treats the viscous term implicitly while the nonlinear advection term 
explicitly is long time stable provided that the time step is sufficiently small in the case with periodic 
boundary conditions. The long time stability in the L 2 and H 1 norms further leads to the convergence 
of the global attractors and invariant measures of the scheme to those of the NSE itself at vanishing 
time step. Both semi-discrete in time and fully discrete schemes with either Galerkin Fourier spectral 
or collocation Fourier spectral methods are considered. 
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1. Introduction. The celebrated Navier-Stokes system for homogeneous incom- 
pressible Newtonian fluids in the vorticity-streamfunction formulation in two dimen- 
sions takes the form 

-Aip = w, 

where uj denotes the vorticity, ip is the streamfunction, and / represents (given) 
external forcing. For simplicity we will assume periodic boundary condition, i.e., the 
domain is a two dimensional torus T 2 , and that all functions have mean zero over the 
torus. 

It is well-known that two dimensional incompressible flows could be extremely 
complicated with possible chaos and turbulent behavior p~3l [HJ EH [27l [38] • Al- 
though some of the features of this turbulent or chaotic behavior may be deduced 
via analytic means, it is widely believed that numerical methods are indispensable 
for obtaining a better understanding of these complicated phenomena. For analytic 
forcing, it is known that the solution is analytic in space (in fact Gevrey class regular 
[12]). and hence Fourier spectral is the obvious choice for spatial discretization. As for 
time discretization, one of the popular schemes [21 [30] is the following semi-implicit 
algorithm, which treats the viscous term implicitly and the nonlinear advection term 
explicitly 



+ V • Vw" - vAuj n+1 = /". (1.2) 
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ere At is the time step, and ui n ,Lo n+1 are the approximations of the vorticity at the 
discrete times nAt, (n + l)At, respectively. The convergence of this scheme on any 
fixed time interval is standard and well-known [TBJ (T7] [T5J [TH] [31] . There are many 
off-the-shelf efficient solvers of (|1.2[) . since it essentially reduces to a Poisson solver at 
each time step. 

It is also well-known that the NSE (jl.lj) is long time enstrophy stable in the sense 
that the enstrophy (i||oj||^ 2 ) is bounded uniformly in time, and it possesses a global 
attractor A and invariant measures [5] [TT] [35]. In fact, it is the long time dynamics 
characterized by the global attractor and invariant measure that are central to the 
understanding of turbulence. Therefore a natural question is if numerical schemes 
such as (|1.2[) can capture the long time dynamics of the NSE (jl.lj) in the sense of 
convergence of global attractors and invariant measures. To say the least, we would 
require that the scheme inherit the long time stability of the NSE. 

There is a long list of works on time discretization of the NSE and related dissipa- 
tive systems that preserve the dissipativity in various forms [31] [32] [9] [TU] [21] [33] [22] 
139] [40] . It has also been discovered recently that if the dissipativity of a dissipative 
system is preserved appropriately, then the numerical scheme would be able to capture 
the long time statistical property of the underlying dissipative system asymptotically, 
in the sense that the invariant measures of the scheme would converge to those of 
the continuous-in-time system |44j . The main purpose of this manuscript is to show 
that the classical scheme (|1.2[) is long time stable in L 2 and H , and that the global 
attractor as well as the invariant measures of the scheme, converge to those of the 
NSE at vanishing time step. 

2. Long time behavior of the semi-discrete scheme. We first recall the 
well-known periodic Sobolev spaces on fi = (0, 2tt) x (0, 27r) with average zero: 

I 4> — and 4> is 27r-periodic in each direction > . (2.1) 
n J 

Hp™ is defined as the dual space of H™ r with the duality induced by the 1? inner 
product. The adoption of £T™ r is well-known [5] [37] since this space is invariant under 
the Navier-Stokes dynamics (jl.lj) , provided that the initial data and the forcing term 
belong to the same space. 

2.1. Long time stability of the scheme. We first prove that the scheme (jl .2[) 
is stable for all time. 

Lemma 2.1. The scheme (|1.2j) forms a dynamical system on L 2 . 

Proof. It is easy to see that for u n 6 L 2 , we have ip n £ H 2 er . Hence V _L ?/' n • 
Wu! n £ Hp^T a for all a £ (0, 1). Therefore, the classical scheme (|1.2j) . which can be 
viewed as a Poisson type problem ui n+1 /At — vAuj n+1 = f — V^ip n ■ Vcj™ + uj n / At £ 
^per a i possesses a unique solution in L 2 (in fact in Hp er ) and the solution depends 
continuously on the data. Therefore it defines a (discrete) semi-group on L 2 . D 

Now we derive the long time stability of the scheme (|1.2j) both in L 2 and in H 1 . 
Our proof relics on a Wente type estimate on the nonlinear term (see Appendix E| . 
which may be of independent interest. 

We first show that the scheme (|1.2[) is uniformly bounded in L 2 , provided that 
the time step is sufficiently small. To this end, we take the scalar product of (|1.2|) 



:=\<P€H m ({l) 
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with 2Atuj n+1 and using the relation 

2(<p-i/;,(p) La = \\(pg-W\\l+\\<p-m (2-2) 



(2-3) 



where || • || 2 denotes the L 2 norm, we obtain 

IK +1 ||2 - IKII2 + \\uj n+1 - uj n \\l + 2 V At\\uj n+1 \\ 2 Hl + 2Atb(^ n ,u: n ,uj n+1 

= 2Ai(r, w " +1 ) L2 

where 

b(il>,u,S>) := (V X V ' Vw,£) i2 = -b{i/),Q,u), (2.4) 

the last equality obtaining upon integration by parts. Using the Cauchy-Schwarz and 
the Poincare inequalities, we majorize the right-hand side of (|2.3p by 

2At||/"|| 2 || w " +1 || 2 < 2c At\\f n \\ 2 \\u J n + 1 \\ Hl < vAt\\u n+l \\ 2 Hl +^At\\f n \\l (2.5) 

Using the Wente type estimate (|A.2j) . we bound the nonlinear term as 
2Atb(ifj n ,uj n ,uj n+1 ) = 2Atb(i/j n ,uj n+1 ,u n+1 -u") 

<2c w At\\v x r\\ H 4" n+1 y4" n -" n+ % (26) 

< I|| w "+ 1 - w "||2 + 2C2Ai 2 ||V ± ^ l ||^ 1 || W "+ 1 ||^ 1 

< - u n \\l + 2C 2 w At 2 \\^f 2 \\uj n + 1 \\ 2 Hl . 

Relations $L3§- ([221) imply 



j n+1 \\ 2 2 - \\u n \\ 2 + - u n \\l + (y- 2ClAt\\Lu n \\ 2 2 )At ||w™ +1 |' 2 



H 1 

r- (2-7) 



At||/ ., 
V 

Here and in what follows, C and c denote generic constants whose value may not be 
the same each time they appear. Numbered constants, e.g., c 42 , have fixed values. 
We are now able to prove the following: 

Lemma 2.2. Let ujq e L 2 and let oj n be the solution of the numerical scheme 
(EH!) . Also, let f G L°°(R + ;iJ) and set H/H^ := ||/||l~(r +; h)- ^Tien t/iere exists 
M = M (||o;o||2,i / , ll/Hoo) such that if 

M< -^ (2 - 8) 

then 

IK|| 2 <M , Vn>0, (2.9) 



< 1 



2r 2 



At 



IKI 



2r 4 

•31/11 



, 1 + S! At , 



,Vn > 0, (2.10) 
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and 

m 2 

^AtglKII^^H^IIl + ^ll/II^Cm-i + lJAt, Vt = l,-..,m. (2.11) 



Proof. We will first prove (|2.10|) by induction on n. It is clear that (|2.10|) holds 
for n = 0. Assuming that (12.101) holds for n = 0, • • • , m, we then have that (|2.9[) holds 
for n = 0, • • • , m, where 

M 2 = Af 2 (|| Wo || 2 , i/, H/IU) - Ml + ^fll/IlL (2.12) 



Then dUTJ) and yield 



l I I 17X1 I 1 9 ^ A . I I « _L 1 1 1 9 . Q 



2 



2" " z 2 " IIJ1 - v 



M\\uj n+1 \\ 2 m < -Atlirila (2-13) 



for all n = 0, • • • , m. Using again the Poincare inequality, the above inequality implies 

\\uj n+1 \\j< + ^At\\r\\l (2.14) 

a av 

where 



a=l+-^At. (2.15) 



-o 

Using recursively (|2.14l) . we find 



m+l 



OJ 



, '^<^M\\l + ^Y,^\\r +1 - i \\l 



v ' — ' or 

1=1 (2.16) 



)— m— 1 r, 4 

W5 + ^ll/ll« 



1 - ,1+ 2^ At 



-in— l-i 



and thus (|2.10p holds for n = m + 1. We therefore have that (|2.10[) holds for n > 
and (|2.9|) follows right away. 

Now adding inequalities (]2.13[) with n from i to m and dropping some positive 
terms, we find 



n 1 1 2 
2 



(2.17) 



<||^ + ^||/||^(m-i + l)Ai, 

which is exactly (|2. 1 1|) . This completes the proof of Lemma 12.21 □ 
Corollary 1. If 

0<At< mini— ^1 =: fc , (2.18) 



then 



\4C72M 2 ' !/ 



<2p 2 , VnAi >T (||wo|| 2j H/lk) := ^ In fJNk ) , (2.19) 

v \ Po 
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where p := (V^cg/^H/Hoo- 

Proof. From the bound (|2.10[) on ||w n |||, we infer that 

IKII|< (l+^At)" n ||u;o||^+^, 

and using assumption (|2.18p on At and the fact that 1 + x > exp(x/2) if a; € (0, 1), 
we obtain 



.nil 2 



For nAt > To, the last inequality implies the conclusion (|2.19|) of the Corollary. □ 

Now we show that the H 1 norm is also bounded uniformly in time under the 
same kind of constraint as for the L 2 estimate. To this end, we first prove that uf 1 is 
bounded for n < N, for some N, and then, with the aid of a version of the discrete 
uniform Gronwall lemma, we show that u> n is bounded for all n > N. 

More precisely, we have the following: 

Lemma 2.3. Let luq G L 2 and let ui n be the solution of the numerical scheme 
Also, let At < ko, with ko as in Corollary{Jl and let r > Scq/u be arbitrarily 
fixed. Then, for n = 1, • • • , No + N r — 1, 

||«l| r t<4(^/^^)(|| Wb ||i P + ^ll/l^) (2.20) 

where Nq = |_7o/At|, with N r — [r/At\ and T that in Corollary^ 

Proof. Taking the scalar product of (|1.2I) with — 2 At Auj n+1 , we obtain 



\\" n+1 \\m - \\" n \\m + - u n \\m + 2vAt\\Aaj n+1 \\ 2 

~2Atb{^ n ,uj n ,Auj n+1 ) = -2At(/",Ao;" +1 ) L 2. 

We bound the right-hand side of (|2 . 2 1 1) using the Cauchy-Schwarz inequality, 



(2.21) 



-2Ai(/«, Au/ l+1 ) L 2 < 2Ai||/"|| 2 ||Au/ l+1 || 2 < ^At\\Auj n+1 \\ 2 2 + ^-At\\f n \\ 2 2 . (2.22) 

Using the Wente type estimate (IA.2|) . we bound the nonlinear term as 

2At&(V>",u/",Aw" +1 ) = 2At6(V> n ,u/ l -w n+1 ,Aw" +1 ) 

+ 2Atb(ip n ,cj n+1 ,Aw n+1 ) 
<2C w At\\V^i J n \\ m \\uj n + 1 - u n \\ m \\ALo n+1 \\ 2 
+ 2C w At\\V^^ n \\ m \\u: n+1 \\ m \\ALo n+1 \\ 2 

< \\\^ n+1 -^ n \\m +2ClAt 2 \\vH n \\ 2 H A\^ n+1 \\l 

v 9C 2 
+ -At\\Au n+1 \\ 2 2 + -^At||vV l |lHi|K +1 ||^i 

< \\\" n+l ^H 2 ffl + 2C? At 2 \\u n \\ 2 2 \\ Au n + l \\l 

v 9C 2 
+ -At\\ALu n+1 \\ 2 + -^Ai||w ,i ||^||w' i+1 ||^ 1 . (2.23) 
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Relations (|!OTj) - (|%23"l) imply 



2C. 



^\\u> n \tiAt\\\u 



n+li|2 



,«+! _ ,.,™ll 2 



IH 1 



- 2C 2 W AtM 2 ) At||A w " +1 || 2 < -Atlirill, 



(2.24) 



from which we find 



where 



a = 1 



2C, 2 



AtM£ > 0. 



(2.25) 



(2.26) 



Using recursively (|2.25p . we find 

1 9 " +1 1 

IK +1 || 2 , 1 <^ T ||c„|| 2 , 1 + -At|| / || 2 <) ^i- 

(X * v — n ' 

2a 2 



< 1 



-AtMn 



ClMl 



l/l 



(2.27) 



Since 2C 2 M^Ai/f < 1/2 by hypothesis (f2~18|) and 

1 - a: > A~ x if x € (0,1/2), 



relation (|2.27[) gives conclusion (|2.20|) of Lemma l2~3l Thus, the lemma is proved. □ 

In order to obtain a uniform bound valid for n > iVo + N r , we need the following 
discrete uniform Gronwall lemma, which has been proved in |39j and we repeat here 
for convenience. 

Lemma 2.4. We are given At > 0, positive integers n ,ni, and positive sequences 

£n, Vn, Cn Such that 



Atri n+ i < |, Vn > n , 

(1 - A<77, i+ i)£„ + i < f n + AiCn+i, Vn > n . 
Assume also that 

712+711 +1 

Ai ^2 Vn < 0-1, 

71 = 712 
712+711 +1 

At < a 2 , 

71 = 712 
712+711 +1 

At ^ < 0,3, 

71 = 712 

for all U2 > riQ. We then have, 

in+i<(-^- + a 2 )e 4a \ Vn>no+m. 



(2.28) 
(2.29) 



(2.30) 



(2.31) 
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Proof. Let m\ and m 2 be such that n < mj < m 2 < mi + Hi- Using recursively 
fl22HD, we derive 



mitniti ^ mitniti mitniti ^ 

U+n 1+ X< II ir^-U-l+A* 2 C« II T3A^- ( 2 - 32 ) 



Using the fact that 1 — x > e 4x , Vx e (0, |) , and recalling assumptions (|2.28[) . and 
the first and second conditions in (|2.30|) . we obtain 

Cmi+ni+l < (£ro 2 -l +fl2)e 4ai . 

Multiplying this inequality by At, summing m 2 from mj to mi + n\ and using the 
third assumption (|2.30l) gives conclusion (I2.31[) of the lemma. □ 

We are now able to derive a uniform bound for ||u; n ||#i valid for sufficiently large 
n. More precisely, we have the following: 

Lemma 2.5. Let luq € L 2 and let uj n be the solution of the numerical scheme 
(HHP- Also, let At < fco, with ko as in Corollary [7J Then there exist constants 
Mi = Mx(u, H/lloo) aradJV = iV(||wo||2,i',||/||oo) such that 

\\uj n \\ H i < Mi, Vn > N. (2.33) 



Proof. Let At be as in the hypothesis, To be as in Corollary [T] r as in Lemma 
Hand set iV := [T a /At\. We will apply LemmaElto (|2~24|) . with £„ = ||w"||^i, 
Tfe = 2C.2 llw"- 1 ^/^, C„ = 2||/||^>, n Q = N Q + 2, m = N r - 2. For n 2 > n Q , we 

12 
2 



compute (taking into account that, by (|2.19p . < 2/5q, for n > No) 



"2+"l + l 9r , 2 



V 



"2 +111 + 1 

At 51 ^ = A< Z 

n— ri2 n—n^ 
n 2 +ni + l na+ni+l „ ~ 

At £ c„ = At £ -ii/ii^-ii/ii: 

^ ' t 7/ 7/ 



J. 



-p r := ai, 



T := a 2 , 



n=n 2 

"2+11 + 1 



n=Tt2 
Ti2+ni +1 



At Cn=At J2 



\ 2 m (by(EH|) 



(2.34) 
(2.35) 
(2.36) 



<- l\\w 

v 



ra 2 -l||2 , H) [ 



||/||^(m+2)At (byQ) (2-37) 



2 

< - 

v 



a 3 . 



By (j23Tj) . we obtain 



4 f2p 2 o , 1 



- ^ + i-H/ll^ +-\\f\\Lr 
v \ r v\\ 



exp 



/16C 2 



Taking N = Nq + N r , we obtain conclusion (|2.33p of Lemma l2"T5l □ 



(2.38) 

(2.39) 
(2.40) 



8 



S. GOTTLIEB ET AL. 



We can summarize the above results in the following: 

Theorem 2.6. The classical scheme il.ty) defines a discrete dynamical system 
on L 2 that is long time stable in both L 2 and H 1 norms. More precisely, for any 
ojq G L 2 , there exist constants k = fco(||wo||2, v, ||/||oo), M Q = A/ n ( 1 1 <^o 1 1 2 , ll/lloo), 
Mi =M 1 (v,\\f\\ 00 ) ondJV = JV(||wo||2,i/,||/||oo) such that 

||w"||a < -Mb, Vn > 0,Vfc e (0,fc ), (2.41) 
\\u n \\ H i < M x , Vn > N,Vk e (0,fc o ). (2.42) 



2.2. Convergence of long time statistics. Here we show that, with time- 
independent forcing, the long time statistical properties as well as the global attractors 
of the scheme (|1.2|) converge to that of the Navier-Stokes sytem at vanishing 

time step size. This is a straightforward application of the abstract convergence result 
(Prop. [2]) in Appendix IB] which itself is a slight modification of the results presented 
in gl]. 

Theorem 2.7. Let d f f = 0. The global attractor and the long time statistical 
properties of the classical scheme U.ty) converge to that of the Navier-Stokes system 
hl-1}) at vanishing time step. 

Proof. We use the abstract convergence result Prop. [2 taking X = B(0, H/l^/^), 
i.e. a ball in L 2 centered at the origin with radius H/l^/^- (The size of the ball needs 
to be adjusted depending on the absorbing property of the scheme.) 

The uniform continuity (H5) of the Navier-Stokes system (jl.ip is a classical result 
[5J |37]. The uniform dissipativity (H3) of the scheme (|1.2p for small enough time 
step with the choice of the phase space X follows from Theorem 12.61 The uniform 
convergence on finite time interval (H4) is proved in Lemma 12.81 below. □ 

Lemma 2.8. Let u be the solution of the continuous system (II. 1| with uj(0) = 
loq G A and uj n that of (|1.2[) with lu° = loq. Assume that f is sufficiently smooth so 
that 

My ~ sup (IIMl-! + IMIIallfitwHia) < 00, (2.43) 

and that Theorem 1,2. 61 holds. Then for At < ko one has 

Ho/ 1 -u(nAt)\\l < AtC(M ,M v ;v) (2.44) 
for all < nAt < 1. 

Proof. We follow the approach in [28, §17] and take d t f = 0. For notational 
convenience, we write t n := nAt and w ra := oj(nAt). Using the identity 

/.(n+l)At 

/ (t - nAt) d tt uj(t) dt = Atd t u\ (n+1)At -Wn+i+Wn, (2.45) 

J nAt 

we have 

Wn+ A~t ~ + V±qJjn ' Vw " ~ " AtJn+1 = f + Rn+1 ■ ( 2 ' 46 ) 
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Here —Aij) n :— u n and the local truncation error is 

-Rn+i := V- L (5?A„+rVw„-V- L Vn+rV<5w n+ i + — / (t-nAt) d tt uj(t) dt (2.47) 

m JnAt 



with 



/-(n+l)At 

Suj n+1 := w n+1 -w„ = / d t uj(t) dt and - A^n+i := ^n+i- (2.48) 

JnAt 

We now consider the error e™ := ui n — uj n , which satisfies 



= -V X V» • Ve" - V^™ • Vw" + i?„+i 

with e° = and -A0" := e". Multiplying by 2Aie™ +1 , we find 

||e" +1 ||2 - ||e«||2 + || e " +1 - e"|ll + 2^|| e "+ 1 || 2 /1 

+ 2Atb(^ n ,e n+ \e n+1 - e") + 2Ai 6(0" , u/\ e" +1 ) (2.50) 
= 2At(i? Il+l! e n+1 ). 

Bounding the nonlinear terms as 

2At (V^n • Ve n+1 , e" +1 - e") < ||e n+1 - e n \\j + At 2 \\ V^V™ ■ Ve rl+1 ||l 

<||e"+ 1 -ell + CSAt 2 K||l||Ve"+ 1 ||l ^ 

where (|A.2|) has been used for the second inequality, and 



2A*(V-V i • Vw",e n+1 ) < 2Ai IjV^" • Ve" +1 || 2 || 
<2C^Ai|| e "|| 2 ||e" +1 || Hl | 



f 12 At 

<uAt\\e n+1 \\^ + ^^\\u> n \\l\\e n f 2 , 



(2.52) 



v 

we obtain, noting that At < ko implies v — C^At \\unW2 > v/2 > 0, 

||e" +1 || 2 + A< (^-C 2 Ai|K|| 2 )|| e " +1 || 2 ffl 

/ r 2 A^ \ (2.53) 
< fl + ^E|| w *|||)|| e n||| + ^11^1^. 



It remains to bound R n +i in H , so for the second term in (|2.47[) we compute, 
for any fixed ip € H , 



\b(ip n+1 ,d t u,ip) \ = \(y ± f ■ Vil) n+ \,d t u) 



L- 



< IkllifillV'n+lll^ll^wll 



(2.54) 



L- 



where (| A.2[) and the identity b(p,q,r) — b(q,r,p) = b(r,p,q) have been used. Similarly, 
for the first term, 



\b(ui n ,d t tp,ip) \ = |(V"V • Vd t ip,uj n ) 



L 2 



(2.55) 

< C w \\v\\m\\Un\\ L A\ d tll>\\H*- 
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The last term in (|2.47[) is readily bounded, and we have by Cauchy-Schwarz, 

(2.56) 



iRn+ifn-i <cAt sup \Ht)\\h \\dMt)f L 2 dt 

te[nAt,(n+l)At] J nAt 



/■(n+l)At 

At / ||flttw(t)H!r-i dt. 

JnAt 



The following bound then follows easily 



/ cAt \ n+1 v-^ 

■!!■]= ||e n+1 ||^ <c(l + — M 2 ) A *H^+illff- 



,n+l 1 1 2 



3=0 



< cAt 2 exp( C ( " + 1)At M'q ) M 2 ((n + l)At) (2.57) 



where 



t ft 

'M|2 



M 2 (t):= / ||9«u;(t0lllr- 1 dt / + sup HOIIi* / I^MOU!* < (2.58) 
Jo t'e[o,t] Jo 

and with it the lemma. □ 

3. Galerkin Fourier spectral approximation. This section is devoted to the 
long time stability of the following Galerkin Fourier spectral approximation of the two 
dimensional Navier-Stokes equations 

71-1-1 n 

UJN At UJN + Pn^Hn ■ Va&) - vAu n N +1 = P N {P). (3.1) 

where w]^, ip^ £ Vn ■— {ah trigonometric functions on Q with frequency in each 
direction at most N}. Pjv is defined as the orthogonal projection from L 2 (ft) onto 
V N . 

Just like for the semi-discrete scheme (|1.2j) . we can show that the scheme (|3.1[) 
is uniformly bounded in L 2 , provided that the time step is sufficiently small. More 
precisely, we have the following: 

Lemma 3.1. Let lu € L 2 and let uj^ be the solution of the numerical scheme 
\3.1\) . Also, let f G L oc (H + ; H) and set \\f\\oo '■— \\f\\L°°(R + ;H)- Then there exists 
M = M (\\uj Q \\2,v,\\f\\oo) such that if 

then 

M|2<M 0! Vn>0, (3.3) 



IKIIa< 1 + 



24 



At 



1 1 1 1 1 - 



/II 



1-1 



2r 2 



At 



, Vn> 0, (3.4) 



and 



7 AMI// 1 



<I|W 



AT 



^II/IIL(™ 



i + l)Ai, Vi = l, 



(3.5) 
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Proof. Taking the scalar product of (|3.ip with 2 At oj]^ we obtain 



n+H|2 ||,..nna , ll,_,n+l_ w n||2 +2l/At || w „+l||| i 



V JV II 2 



+ 2Ai = 2At(/», u^+V . 



(3.6) 



Using the Cauchy-Schwarz inequality and the Poincare inequality, we have the fol- 
lowing bound for the right-hand side of (|3.6j) : 



2At(r,^ +1 ) L 2 < 2At||/»|| 2 |K +1 || 2 < 2Atc \\f n \\ 2 \\uj^ +1 \\ m 



pn\\2 
2j 



(3.7) 



< vAt\\wl +l f H1 + -fAtWf 

whereas the nonlinear term can be bounded using the Wente type inequality (|A.2[) as 

2At b(^,cj%,LJ n N +1 ) = 2At b{r N ^ n +\oj n + l - u n N ) 



<2C w At\\V x ^\\^\\w n N +1 \\ H i\\w 



N 



^ I n+1 n II 2 



2 
1, 



N 

n+1 
N 



, vll2T 2^Ai 2 ||V J -^ll^ll^ +1 |l^ 



(3.8) 



J 7\' 1 1 2 



2ClAt 2 \\co 



n ||2||, ,n+l||2 



N\\2\\" J N 



H 1 ■ 



Relations (|3l|-(|3JiJ) imply 



n 1 1 2 
Nl\2 



,n+l 



£ + (i/-2C£Af|K||^)Af| 



,n + l |2 



< -Ai||/ n || 2 . 

V 

By induction, one can prove that if At satisfies (|3.2j) . then 



(3.9) 



\\^\\l< I 1 



< 1 



2c 2 



2c 2 



< IMI! 



At 
At 

24 



2c, 



J N\\2 



i ii/n 



i- i 



INI 



2q 



2c 2 



At 



2cl 



At 



11/11 



=: A/ 2 (|| Wo ||2,M/||oo), Vn>0, 



(3.10) 



from which conclusions 
Adding inequalities 
the time restriction 



and (I3.4j) of the Lemma follow right away, 
with n from i to to and recalling the bound 
wc find 



and 



-Ai^Hi. 

n— z 



< IKI 



i II 2 
2 



^At 



11/ 



n||2 
2 



(3.11) 



<IKvll! + -ll/llL(™-* + i)Ai, 



(3.12) 



which is exactly conclusion (I3.5[) of Lemma 13.11 This completes the proof of the 
lemma. □ 
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From bound (|3.4p we can also derive the following 
Corollary 2. If 

0<At< mini— * M\ =: fc , (3.13) 



\4C 2 M 2 ' * 



then 



K||l<2pg, VnAt>T (|M|aJ|/||oo):=^ Infill, CS.N-i 



w/iere p := 11/11 



Using Lemma 12.41 we can prove a result similar to Lemma 12.51 More precisely, 
we have the following: 

Lemma 3.2. Let ojq G L 2 and let lu^ be the solution of the numerical scheme 
(3.1]) . Also, let At < fco, with kg as in Corollary [Jl Then there exist constants 
Mi =M 1 (!/,||/|| 0o ),iV = AT(||wo||2,^||/||oo) such that 

\\unWh 1 < M u Vn > N. (3.15) 



Proof. Taking the scalar product of ([57T]l with -2AtAu;™ +1 , we obtain 

- 2Atb(^,w%, Aw n N +1 ) = -2At(f n , Auj n N +1 ) L 2 . (3.16) 
Using the Cauchy-Schwarz inequality, we bound the right-hand side of (|3.16[) as 

-2Ai(/ n , Au n N +1 ) L . < 2At||r|| 2 ||A^ +1 || 2 < ^At\\ALj% +1 \\l+*At\\f n \\l (3.17) 

Using the Wente type estimate (IA.2I) . the nonlinear term can be bounded as 

2Ai6(V^,< ) Au™+ 1 ) = 2Ai&(^ 

<2C w At\\V ± ^\\ m \\^ +1 - ^\\ m \\A^ +1 \\ 2 



-2C tu Ai||V ± V^lk 1 ||^ +1 |k 1 ||A^+ 1 || 2 



< + 2C 2 Ai 2 ||V^|| 2 , 1 ||Ac^ +1 || 2 



v . . „ . „ 1 1 ,,o 2C 2 



+~At\\A^\\l + ^At||V>£||^|K +1 ||^ 
1, 



< - w ft||^ + 2CiA^||^|| 2 || Ac^+12 



^ A . i I a I j -!- I I I '"> ■ ^ 



+-At||A w ™ +1 ||| + ^Ai||^|||||^ +1 || 2 , 1 . (3.18) 



Relations (l3J6)) - (j3A8)) imply 



or 2 \ 1 

—m.]]u^r\\lAt ll^llix - llc^ll 2 , + -Hw^-w^ll? 



2C 2 AiM 2 ) At||Ac^ +1 || 2 < ^Atlini 2 , 



(3.19) 
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from which we find 



where 



^Ulr^ilKII^ + ^-Atll/ll 



9C 2 

a=l- -^AtM* > 0. 



(3.20) 



(3.21) 



Now let No = [To/ At], with To being given in Corollary [21 and for r > 8cq/v 
arbitrarily fixed, let N r = [r/At\. We are going to apply Lemma \2. 41 to (|3.19j) . with 
e„ = IK^i, r, n = 2Cl IK-'HIM Cn = 211/H^/i/, n = iVo + 2, n x = JV r - 2. For 
n>2 > no, we compute (taking into account that, by (|3 . 14[) . Hw^Hf < 2pp, for n > No): 



(3.22) 
(3.23) 
(3.24) 



n 2 +rai+l 


n 2 +ni+l 


9^2 

W 10 || ,71 — 


AC 2 


At £ ^ 


= At £ 




n=«2 


n2+ni+l 


H2+ni+l 






At £ Cn 


= At £ 


-llfll 2 < 
1/ - 


-JfWl.r :=a 2 


n=n2 


n=n 2 


"2+ni+l 


n2 +m + l 






At £ e„ 


= At £ 




(by (33}) 


n— n2 


71 — 712 







< - I llw 

V 

2 

< - 

V 



,n 2 -l 1 1 2 _,_ L 



5 + -H/IISo(ni+2)At (by (EH) (3.25) 



2 Po + ^ 



L n ii ...2 



By (|2.3ip . we obtain 



J N\\m 



< 



=: a 3 . 



+ -H/ll»r 



exp 



16C3 



Taking N = N 



= : Af?(i/, H/Hoo), Vn>iV + AT r . 
N r , we obtain conclusion (|3. 15|) of Lemma [3T2l □ 



(3.26) 

(3.27) 
(3.28) 



4. Collocation Fourier spectral approximation. Here we consider the collo- 
cation Fourier spectral spatial approximation of the scheme (|1.2j) . In order to maintain 
the long time stability of the fully discretized scheme, a common technique of using 
a modified form of the nonlinear term is utilized (see for instance [3S]). Moreover, 
we will use an alternative approach for the nonlinear analysis: instead of apply- 
ing the Wente type estimate, we will use ||V^||l°o, which is in turn bounded by 
llV'llj^allV'll/p 6 ! G (0, 1). This alternative approach leads to a slightly more restric- 
tive time step restriction for stability, but has the advantage of easy adaptance to the 
fully discrete collocation Fourier approximation. 

4.1. Fourier collocation spectral differentiation. Consider a 2-D domain 
f2 = (0, L x ) X (0, L y ). For simplicity of presentation we assume that L x = L y = Lo = 1 
and L x = N x -h Xl L y — N y ■ h y for some mesh sizes h x = h y = h > and some positive 
integers N x = N y = 2N + 1. All variables are evaluated at the regular numerical grid 
(xi,yj), with Xi = ih, yj = jh, < i,j < N. 
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For a periodic function / over the given 2-D numerical grid, assume its discrete 
Fourier expansion is given by 

[N/2] 

fij= E (fc)jn,h^ i(klXi+llVi) ' (4.1) 

k u h = -[N/2] 

Note that may not be the regular Fourier coefficients, due to the aliasing error. 
In turn, its collocation interpolation operator becomes 

N 

X N f(x)= ]T (fj) kuh e 2 ^ x+l ^. (4.2) 

k!,h=-N 

As a result, its collocation Fourier spectral approximations to first and second order 
partial derivatives (in x direction) are given by 

N 

{V Nx f) itj = J2 {^i){f^) hlM ^ k ^ +l ^\ (4.3) 

k lt li=-N 
[N/2] 

(P^J)ij= E (~^ 2 k 2 )f klM e 2 ^+ l ^l (4.4) 

fci,ii=-[iV/2] 

The corresponding collocation spectral differentiations in y directions can be defined 
in the same way. In turn, the discrete Laplacian, gradient and divergence can be 
denoted as 

Ajv/ = (V 2 Nx + V 2 Ny ) f, \7 N f = ( ® Nxf ), Vjv • ( f f ) = 2Wi + V 



' N 

\ ^NyJ J \ J2 / 

at the point-wise level. 

Moreover, given any periodic grid functions / and g (over the 2-D numerical grid), 
the spectral approximations to the L 2 inner product and L 2 norm are introduced as 

27V 

\\f\\ 2 = VifJ), with (f,g) = h 2 kj9i,r (4-6) 

i,j=a 

Meanwhile, such a discrete L 2 inner product can also be viewed in the Fourier space 
other than in physical space, with the help of Parseval equality: 

N N 



</,<?} = E (fc)k^)^ h = E (AulAfFhiM, (4-7) 
k 1 ,l 1 =-N ki,h=-N 

in which (fc)ki,h , (9c)ki,h are the Fourier interpolation coefficients of the grid func- 
tions / and g in the expansion as in (14.1[) . Furthermore, a detailed calculation shows 
that the following formulas of summation by parts are also valid at the discrete level: 

(f, V ^ ' ( ? )) = - (^Nf, ( I 1 ) ) , (/, A N g) = - (V N f, V N g) . (4.8) 
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4.1.1. A preliminary estimate in Fourier collocation spectral space. It 

is well-known that the existence of aliasing error in the nonlinear term poses a serious 
challenge in the numerical analysis of Fourier collocation spectral scheme. To over- 
come a key difficulty associated with the H m bound of the nonlinear term obtained 
by collocation interpolation, the following lemma is introduced. The result is cited 
from a recent work [15) . and the detailed proof is skipped. 
Lemma 4.1. For any <p e P 2N in dimension d, we have 

\\lN^\\ Hk <(V2) d M Hk . (4.9) 



In fact, an estimate for the k = case was reported in E's work [7J [8], with 

the constant given by 3 d , while this lemma sharpens the constant to %/2 . The 
case with k > | = 1 was covered in a classical approximation estimate for spectral 
expansions and interpolations in Sobolev spaces, reported by Canuto and Quarteroni 
[3J. However, due to the additional regularity requirement for interpolation operator 
analysis, the case of k = 1 was not covered in any existing literature, which we require 
for the H 1 bound of the nonlinear expansion in the global in time analysis. 

4.2. The first order semi-implicit scheme. The fully discrete pseudo-spectral 
scheme follows the semi-implicit idea of (jl.2p and (|3.1[) : 



— + - (u".VW + V w • (u n u n )) = uA N u, n+1 + f\ (4.10) 

-A N ip n+1 =u n+1 , (4.11) 
u n+1 = VU n+1 = {V Ny r + \V Nx r +1 ) ■ (4.12) 

It is observed that the numerical velocity u n+1 — Vj^'ip n+1 is automatically 
divergence- free : 

Vat • u = V Nx u + T> Ny v = -T> Nx (V N yip) + T> Ny (T> Nx ip) = 0, (4.13) 

at any time step. Meanwhile, note that the nonlinear term is a spectral approximation 
to |m"-Vw and • (uw) at time step t n . Furthermore, a careful application of 
summation by parts formula (14.81) gives 

(u, u-\7 n uj + Vat • (uw)) = (u, u- Vatw) — (Vjvw, uui) = 0. (4-14) 

In other words, the nonlinear convection term appearing in the numerical scheme 
(|4.10[) . so-called skew symmetric form, makes the nonlinear term orthogonal to the 
vorticity field in the L? space, without considering the temporal discretization. This 
property is crucial in the stability analysis for the Fourier collocation spectral scheme 

In addition, we denote U n = (U n ,V n ), u> n and ip n as the continuous versions 
of u n , u n and ip n , respectively, with the formula given by (|4.2p . It is clear that 
U n ,u n ,ij> n £ P N and the kinematic equation Ai/>™ = oj n , U n = V^ip n is satisfied 
at the continuous level. Because of these kinematic equations, an application of elliptic 
regularity shows that 

W\\ Hm+ , <C\\^\\ Hm , ||V"|| Hm+2+Q < C\\u, n \\ Hm+a , (4.15) 
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in which we used the fact that all profiles have mean zero over the domain: 



V>™ = 0, U n = [-d y tp r \ d x ip n ) = 0, = AV>" = 0. (4.16) 

Moreover, it is clear that the Poincare inequality and elliptic regularity can be applied 
because of this property. 

Lemma 4.2. Let luq £ L 2 and let ui n be the solution of the numerical scheme 
UAd\) - UJ^) . Also, let f £ L°°(B. + ;H) and set H/H^ := ||/||ioo (R+;H) . Then there 
exists Mo = Mo(||wo||2, ||/||oo) such that if 

U -W ( " 7) 

then 

< M , Vn > 0, (4.18) 



2 



< 1 



2c 2 



At 



2 



9r A 

— Ilfll 2 

v i ii J Hoc 



1+ 2? A * 



, Vn > 0, 



(4.19) 



and 



m 2 

faY, \\u, n f H * < WuFfv + ^\\f\\l(rn - i + l)At, Vi = 1, • • ■ ,m.(4.20) 

n—i 



The proof of this lemma is organized as follows. First, an H s a-priori assumption 
for the numerical solution u n is made. In turn, this assumption leads to a global 
in time L 2 bound, with a standard application of Sobolev embedding and Holder's 
inequality. However, this L 2 bound is not sufficient to recover the a-priori assump- 
tion, due to the fact that the Wente type analysis is not available for the collocation 
spectral approximation. Instead, a global in time H 1 stability can also be derived 
with the help of the leading L 2 bound. Moreover, both the global in time L 2 and H 1 
bound constants are independent of the a-priori constant C\. As a result, the a-priori 
assumption can be recovered so that an induction can be applied to established the 
above lemma. 

4.3. Leading estimate: L°°(0, T; L 2 ) n L 2 (0, T; H 1 ) estimate for u. Assume 
a-priori that 

ll w "ll_ff 5 < Ci, u)" is the continuous version of w™, (4-21) 

for some S > at time step t n . Note that C\ is a global constant in time. We are 
going to prove that such a bound for the numerical solution is also available at time 
step t n+1 . 

Taking the discrete inner product of (|4.10p with 2Atu n+1 gives 

\\io n+1 \\ 2 - \\Lu n \\ 2 + \\tu n+1 - w n \\l + 2isAt\\V N Lu n+1 \\ 2 

= -At (u n -V N Lu n + Vat • (u n w n ) ,u/ l+1 ) + 2At (f\uj n+1 ) , (4.22) 
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(4.26) 



in which the summation by parts formula (I4.8[) was applied to the diffusion term. A 
bound for the outer force term is straightforward: 

2(r> n+1 ) < 2 ||/"|| 2 ■ || w " +1 || 2 < 2C7 2 ||r|| 2 • ||V^" +1 || 2 

< P^X + — < ~, P^t + ^,(4-23) 

in which a Poincare inequality 

||^" +1 || 2 <C , 2 ||V A r W n+1 || 2 , (4.24) 

was used in the third step. For the nonlinear term, we start with the following 
rewritten form: 

-At{u n -V N uj n + V N ■ {u n uj n ) ,uj 11+1 ) 
= -At{u n -V N uj n+1 + Vat • (u n uj n+1 ) 

+At (u n -V N {Lu n+1 - u n ) + V N ■ (u n {uj n+1 - io n )) , w n+1 ) . (4.25) 

The first term disappears, using a similar analysis as ()4. 14|) : 

(u' l -V N UJ n+1 + V N ■ (V 1 U/ 1+1 ) ,LU n+1 ) 

= {io n+l 1 u n -\/ N uj n+1 ) - (\7 N uj n+ \u n uj n+1 ) = 0. 
For the second term, the summation by parts formula (|4.8[) can be applied: 

/ ,n v-7 ( n+1 n\ n+l\ / n+1 n t-7 /„.n n+l\\ /a 

(u ■ Vat(w — w )j w > = — ( w — w , Vat • (it w )) , (4.27) 

/w ^« n/ n+1 n\\ n+l\ / .n+1 n „ n v7 n+l\ /,i oo\ 

^Vjv • (M (w — uj )) , w ) = — i 10 — lu ,u ■ VnLli ) , (4-28) 

For the term Vat • (u n u n+1 ), we note that it cannot be expanded as u n ■ Vatw™ +1 , 
as in the Fourier-Galerkin approximation, even though u n is divergence-free at the 
discrete level ()4.13j) . In the collocation space, we have to start from 

Vat ■ (u n w n+1 ) = V Nx (u n oj n+1 ) + V Nv (v n u n+1 ). (4.29) 

To obtain an estimate of these nonlinear expansions, we recall that U n = (U n ,V n ), 
u) n+1 and t/>" +1 are the continuous versions of u n , ui n+1 and ip n+1 , respectively. Since 
U n , u n+1 £ P N , we have U n u) n+1 £ P 2N and an application of Lemma 14.11 indicates 
that 

||Pa^KV" +1 )|| 2 = \\d x l N (U n cj n+1 )\\ 2 < 2||^(C/W l + 1 )|| 2 , 
\\V Ny (v n ui n+1 )\\ 2 = \\d y l N (V n oj n+1 )\\ 2 < 2 \\d y (V n cv n+1 )\\ 2 . (4.30) 

Subsequently, a detailed expansion in the continuous space and an application of 
Holder's inequality show that 

\\d x (u n u} n+1 )\\ 2 = ||J7> ,1+1 + j/"< +1 || 2 < ||t/>" +1 || 2 + \\u n ^ +1 \\ 2 

< \\U2hw-» ■ \\" n+1 \\ L vs + \\U n \\ L ~ ■ \H +1 \\ 2 ■ (4-31) 
Furthermore, a 2-D Sobolev embedding gives 

\\" n+1 \\ L v* <C\\Un H * \\" n+1 \\ m <0\\u: n \\ HS ||Vu;" +1 || 2 (4.32) 
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in which the elliptic regularity (|4.15[) and the Poincare inequality were utilized in the 
last step. The second part in (|4.31l) can be handled in a straightforward way: 

\\U n \\ LOO ■ ||^ +1 || 2 < C\\U n \\ Hl+s ■ \\Vu n+ \ < C\\u: n \\ HS ■ ||Va." +1 || 2 ,(4.33) 

with the the elliptic regularity (|4.15[) applied again in the second step. A combination 
of (l4~52l and yields 

\\d x (U n u> n+1 )\\ 2 <C\\u n \\ H , • ||Va;" +1 || 2 . (4.34) 

Similar estimates can be derived for ||<9j,(U™u;™ +1 ) || 2 . Going back to (|4.30p . we arrive 
at 

||Vjv • {u n Lo n+1 )\\ 2 <C\\u n \\ HS ■ ||Vu;" +1 || 2 =C\\u n \\ HS • ||VW l+1 || 2 , (4.35) 

in which the second step is based on the fact that u> n ,uj n+1 € P N , so that the 
corresponding L 2 and H s norms are equivalent between the continuous projection 
and the discrete version. In addition, the nonlinear term in (I4.28|) can be controlled 
in a similar way: 

||« n • V WW " +1 || 2 < ■ ||V ArW "+ 1 || 2 = C\\u> n \\ HS ■ ||Viva; n+1 || 2 , (4.36) 

with a discrete Sobolev imbedding inequality applied in the second step. Therefore, 
a substitution of (|Q5j| - (jOSll into (14^5) . ||OB)> . QgHh&M results in 

-At{u n -V N uj n + Vjv • {u n uj n ) (4.37) 

< CAt\\u n \\ H t ■ ~ w "|| 2 ' ll^^la 

< CCi&t \\uj n+1 -to n \\ 2 ■ ||VArw' l+1 || 2 



< -vAt\\V N u n+ ^+ < =^At\\u n+l -u n \\l. (4.38) 
Its combination with ([OS) , (EL"2"6"|) and P~2"2"j) leads to 



n+l 1 1 2 II, ,n\\2 



w »||3 + I i _ 9^±At) ||^+ 1 - w «||2 + l ,Ai||V Ar ^+ 1 ||2 < ^i^At 



V 



(4.39) 

Under a constraint for the time step 
G 3 C? . 1 



At < -, i.e., At < =7-, (4.40) 

v "2' " 2C 3 Cl V ; 



we arrive at 



1 



™+l||2 _ II,.,"- II 2 , _||,..n+l _ ,.,ni|2 , „ A . ii . . n+l || 2 



OJ — UJ 



12 



vAt\\^ N uj n+1 \\ 2 2 < C* 4 Ai (4.41) 



with C4 = (2C|M 2 )/j/. Furthermore, an application of the Poincare inequality (|4.24[) 
implies that 



j n+1 \\l - \\u n \\l + C b vAt\\uj n+1 \\ 2 2 < C 4 Ai, with C 5 = i. (4.42) 

^2 
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Applying an induction argument to the above estimate yields 

< (i + evA<)-(" +1 >|| w ||2 + £l 

\\u n+ % < (l + a^At)- (n+1)/2 |k || 2 + y^:-C 6 . (4.43) 

Note that Cq is a time dependent value; however, its time dependence is in exponential 
decay so that a global in time bound is available. 

In addition, we also have the L 2 (0,T;iJ 1 ) bound for the numerical solution: 

N k 

vAt J2 \\ V Nio k \\l < \\u l \\ 2 + d (T* - t) . (4.44) 

k=i+l 

However, it is observed that the a-priori estimate (I4.43|) is not sufficient to bound 
the H s norm (|4.21[) of the vorticity field. In turn, we perform a higher order energy 
estimate L°° (0, T; H 1 ) n L 2 (0, T; H 2 ) for the numerical solution of the vorticity field. 

4.4. L°°{0,t 1 ;H 1 )nL 2 (0,t 1 - 1 H 2 ) estimate for w. Taking the inner product of 
(j4~TTJT) with -2AtA N uj n+1 gives 

||Vjvw n+1 ||2 - \\V N uJ n \\ 2 2 + ||Vjv - u") Hi + 2uAt\\A N io n+1 \\ 2 2 

= At (u n -V N uj n + Vat • (u n uj n ) , A N Lj n+1 ) - 2At (/", A N cj n+1 ) . (4.45) 

The Cauchy inequality can be applied to bound the outer force term: 



-2 (r, aw +i ) < ^ ||a^" +i ||J + 2 - iiriia 

1 ,, . „_ L1 „2 2M 



2 



< l I /||A Jv w n+1 ||; + ^-. (4.46) 

For the nonlinear terms, we first make the following decomposition: 

u n -\7 N uj n = -u n -\7 N (uj n+1 - ui n ) - (u n+1 - u n )-S7 N uj n+1 

+ u n+1 -V N uj n+ \ (4.47) 

Vjv ■ (u n u n ) = Vjv ■ (-u n (uj n+1 - oj n ) - (u n+1 - u n )uj n+1 

+ u n+1 uj n+1 ). (4.48) 

For the first term, the a-priori assumption (|4.21[) gives 

I Tl V7 / n+l ^ ii nil ||t-7 / n+1 n \ 1 1 

||— u -Vat (uj t -w )\\ 2 < \\u ■ || Vjv (w - w J|| 2 



12 

<CCi||V JV (w" +1 -w")|| 2 , (4.49) 



in which we applied the discrete Sobolev inequality in 2-D: ||m ti ||oo < C| |w"|| H i+ 5 < 

h 

C||u; n ||jy«. This in turn leads to 



At(-u n -V N {Lo n+1 - u n ),A N oj n+1 ) 

< CC t At 1 1 Vat (uj n+1 -w n )|| a • HAatu/ 141 

< -^At||A A ,w n+1 ||2 + ^Ai||Viv(w" +1 -w ,l )||2. (4.50) 



2 II Jv ~ 112 

2 
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The conservative nonlinear term Vjy ■ (u n (uj n+1 — w™)) can be analyzed as in (|4.29l) - 

||Vjv ' (u n (uj n+1 - uj n ))\\ 2 < \\D Nx (u n (w n+1 - u n j) |[ 2 + \\V Ny (v n (uj n+1 - w n j) || 2 

< 2 (\\d x (U n (u> n+1 - w"))|| a + \\d y (V n (0J n+1 ~ w n ))|| 2 ) , (4.51) 
\d x (U n {u; n+1 - a;' 1 )) || 2 = \\U™{u> n+1 - u> n ) + £/"(u/ l+1 - w n ) a || 2 

< UuriUa-*) • ||" n+1 -^ILv, + II^Ht- • -«")*l| a 

< C|K|| ffS • ||V(^" +1 -^")|| 2 < CCt \\W N (oj n+1 -uj n )\\ 2 , (4.52) 
\\d y (V n (u> n+1 - u> n )) || 2 < CC 1 ||VArK +1 - w ")|| 2 , (4.53) 

with the help of the elliptic regularity (I4.15[) , Poincare's inequality and 2-D Sobolev 
embedding. Consequently, we see that the first part of the nonlinear term (|4.48[) has 
the same bound as (|4.49|) : 

|| V* • (u n (w n+1 -w n ))|| a < CC X || V w -^")|| 2 , (4.54) 

which in turn leads to an estimate similar to (|4.50|) : 

At (-Vat • (u n (u; n+1 - u n )) , A N uj n+1 ) 

< ^vAt\\A N uj n+1 \\l+ < ^At\\V N (Lj n+1 -u n )\\l. (4.55) 

For the second term in (|4.47|) . we start with the following Sobolev inequality: 



||V^ +1 || 2 - ||Vu," +1 || 2 < ||w" +1 || H1 < C||w n+1 ||2 /2 ■ \\^ n+1 \\f 2 

< C||w n+1 ||2 /2 • ||Aa;" +1 || 2 /2 < CCl' 2 ||Au;™+ 1 || 2 /2 , (4.56) 

in which an elliptic regularity ||o; n+1 1| ff2 < C || Au n+1 1| was utilized in the second 
step and the leading L 2 estimate (|4.43p was used in the last step. Similarly, we also 
observe that the kinematic relationships 

U n+l - U n = V x (i/> n+1 - V'") , A(ip n+1 --0") =cj n+1 -u> n , (4.57) 

indicate the following Sobolev estimates: 

||m" +1 -m"|| < \\U n+1 -U n \\ r 

II 1 1 oo — II 1 1 ^ 

< C \\U n+1 - U n \\ m+S < C \\^ n+1 - ^ n \\ H2+s < C || W » +1 - u n \\ RS 

.n+l , .n II 1-5 II, .ra+1 , .nil 5 

< C \\u> n+1 - u n \\l~ 5 || V (u/ l+1 - w n ) ||* 

< C (2C 6 ) 1_ * || V (u; n+1 - u") ||* , (4.58) 

in which estimate (|4.43[) was used in the last step. Consequently, a combination of 
(|4~56) and (|4~58t indicates that 

||(« n+1 -tt")-VAro; ,l+1 || 2 < \\u n+1 - u* 1 ^ ■ ||Vatu/ 1+1 || 2 

< CCl' 2 (2C & ) 1 - S || V (u," +1 - w")||* • ||Au," +1 || 2 /2 

< CCl 12 (2C 6 )^ S || Vat - io n ) ||* • \\A N u n+1 f 2 ' 2 , (4.59) 
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due to the fact that u> £ P . In turn, the following estimate is obtained 

At {- (u n+1 - u n ) ■V N LJ n+1 , A N uj n+1 ) 

< CCl ,2 At ||Vjv {uj n+1 - Lu n )\\ 5 2 • \\A N uj n+1 f 2 /2 . (4.60) 

Meanwhile, the second conservative nonlinear term in (|4.48[) . Vat-((m™ +1 — u n )uj n+1 ) , 
can be expanded and analyzed in a similar way: 

V/V„. n+l „ n\ n+l\|| ^ ||t-» / / n\ n+l\|| , ||t-» / / n+1 n\ n+l\ll 

, at ' ((" - M )w J|| 2 < |PjVa - u V j|| 2 + IPJV» \\ v ~ V ) w J|| 2 

< 2 (||a x - C/")w n+1 ) || 2 + \\d y ((V n+1 - V n )u) n+1 ) || 2 ) , (4.61) 

\\d x ((U n+1 ~ U n )u) n+l )\\ 2 = \\{U n+l - U n ) x uj n+1 + {U n+1 - U n )o 



J x 1 1 2 

< \\(U n+1 - U n ) x \\ L2/(1 _ S) ■ ||o; n+1 || L2/a + \\U n+1 - U n \\ Loc ■ \\^ +1 \\ 2 

<C\\U n+1 -u n \\ H1+5 • ||v^ Tl+1 || 2 

< CCl 12 p N (cj n+1 - u n )\\ 5 2 ■ ||A ArW n+1 ||J /2 , (4.62) 
\\d y ({V n+l - V n )uj n+1 )\\ 2 < CCf 2 || Vat (oj n+1 - u n )\\l • ||A A rw ,i+1 ||J /2 . (4.63) 

Again, the elliptic regularity (|4.15p . Poincare's inequality and 2-D Sobolev embedding 
were repeatedly used in the analysis. As a result, its combination with (|4.60[) leads to 

At (- (u n+1 - u n ) -V N uj n+1 - Vat • ((u n+1 - u n )u n+1 ) , A N uj n+1 ) 
< C 7 Cl /2 At ||Vjv (w n+1 - w n )||a • ||A A r W " +1 ||J /2 . (4.64) 

We can always choose < 5 < \ , so that an application of Young's inequality (ab < 
^ + ^ with i + i = 1) gives 

p q p q ' 

|| K +1 - u n ) \\l • || &Nv n+1 \\l /2 < c 8 \\v N K +1 - w ») nf + II a^ 1 !) 2 , 

2CVC 6 

with C 8 = ]( 3 -^] . (4.65) 
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Furthermore, since 45 < 2, we can apply Young's inequality to ||Vjv — u/™) ||2 

and obtain 

C S \\W N - w ») ||f < — ^ HVjv - w «) || 2 + Cfc, (4.66) 

c 7 c 6 

in which Cg depends on C6, C7, Cs and 8. As a result, substituting (|4.65[) ~ (14.66[) into 
(|4.60j) gives an estimate for the second nonlinear term: 



At (- (u n+1 - u n )-V N uj n+1 - V N ■ ((u n+1 - u n )io n+1 ) , A N uj n+1 ) 

1 2 1 „ , 11 A ||2 

2 



< At || Vjv - u n ) H2 + || Aw™ +1 H2 + CioAt (4.67) 



with C10 = C 7 C* 6 3/2 C* 9 

The third nonlinear term in (|4.47l) . (I4.48|) can be analyzed in a similar way. We 
first look at u n+1 ■ Wn^" +1 - A bound for ||m" +1 ||oo can be obtained in the same 
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fashion as (I4.58|) : 



K +1 L <C\\U n+1 \\ H1+5 <C\\rp n+1 \\ H2+5 <C\\uj n+1 \\ H6 <C\\u n+1 \\l ^Hw^ 1 !!*, 
< C||u;™ +1 ||2"^ ||Au/ l+1 ||f < CCg" 1 HAato;'^ 1 !!! . (4.68) 
Its combination with (|4.56[) shows that 

At (u n+1 -\7 N Lo n+1 , A N uj n+1 ) < At ||tt" +1 || oo • ||VArw n+1 || 2 • ||Ajva;" +1 ll 



12 ||~«~ 112 



3 + i5 



<CCs /2 At\\A N oj n+1 \\ 2 2 . (4.69) 
This analysis can be applied to the term Vn ■ (u n+1 us n+1 ) in the same way: 

|| Vat • (u n+1 uj n+1 ) || 2 < \\V Nx (u n+1 oj n+1 ) || 2 + \\V Ny (v n+1 uj n+1 ) || 2 

< 2 (\\8 X (U n+1 u> n+1 ) || 2 + \\d v (V n+1 u> n+1 ) || 2 ) , (4.70) 

\\d x (U n+1 u: n+1 )\\ 2 = \\U x l+1 u> n+1 + U n+1 uj n + 1 \\ 2 

<r llrr ri + 1 ll ILi n + 1 ll x llr/"'^ 1 !! 

-\\ u x ||l=/(i-4) ll w || L 2/s + ||^ \\ L oo ■ \\ U X || 2 

<C\\U n+1 \\ Hl+s • ||Va;" +1 || 2 < CCl /2 \\A N uj n+l \\^ , (4.71) 

\\d y (V n+1 w" +1 )|| 2 < CCl' 2 llAjvo;'^ 1 !!^ . (4.72) 
As a result, we arrive at the following estimate: 

At (u n+1 -\7 N iu n+1 + Vat • (u n+1 uj n+1 ) , A N uj n+1 ) 



< C 11 C 6 3/2 Ad|A A r W "+ 1 || 2 2 . (4.73) 



Again, since < 2, we can apply Young's inequality and obtain 

HAA^Uf < — iL^HAW^lg + Ci* (4.74) 

in which Cn depends on C6, Cn and 5. Going back to (|4.73p . we have an estimate 
for the third nonlinear term: 



At (u n+1 ■ \7 N uj n+1 + Va, • (m"+V 1+1 ) , A N ui n+1 ) 

< ±v£t \A N u n+1 || 2 + C 13 At, (4.75) 

with C13 - C 12 CnCe /2 . 

Finally, a combination of ([£35]) - (|Qg]) , (|OD|l . (|Q5|) . (jlTFfl) and (|4T75l) results in 

||V ArW "+ 1 ||2 - IIV^IH +(l-(l+ ^] At) \\V N (co n+1 - u») \\ 2 



^At|| Aj\ra; Tl+1 H2 < ( — + Cm + C13 ) At. (4.76) 



+±i/Af||A J ^»+ 1 ||*< (^f 
Under a constraint similar to (|4.40p and a trivial constraint At < | for the time step: 

1 At< -, At<-, i.e., At<min , (4.77) 

v "2' -4' " V2C 14 C 1 2 4y' V ' 
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we have 

l|VA^ n+1 ||^ - ||Va^"||2 + ~ ||Vjv (cy i+1 - cj n )\\ 2 2 + ±vAt\\A N oj n+1 \\l < Cu>At, 
2M 2 

with C 15 = — + Cm + Ci 3 . (4.78) 
Furthermore, an application of elliptic regularity 

||VArC^ +1 || 2 <Ci 6 ||A A r W "+ 1 || 2 , (4.79) 

implies that 

\\V N u n+1 \\ 2 2 ~\\V N u n \\ 2 2 + C 17 isAt\\V N u n+1 \\ 2 2 < CisAt, with C17 = =L-. (4.80) 

2o 16 

Applying an induction argument to the above estimate yields 

||Vjvw n+1 ||2 < (l + C 17 vAt)-( n+ V\\V N Lo°\\ 2 + fV ' ' 



Ci 7 v 



\\V N u n+1 \\ 2 < (l + C 17 vAt)- R i 1 \\V N uj \\ 2 + J^- :=C 18 . (4.81) 

Again, Cis is a time dependent value; however, its time dependence is in exponential 
decay so that a global in time bound is available. 

In addition, we also have the L 2 (0,T; H 2 ) bound for the numerical solution: 



-vAt \\ A nco% < llViv^Hl + C u (T* - e) . (4.82) 

k=i+l 



4.5. Recovery of the a-priori H 5 assumption (|4.2ip . With the L°° (0, T; L 2 ) 

and L°°(0,T; H 1 ) estimate for the numerical vorticity solution, namely (|4.43[) and 
(|4.8ip . we are able to recover the H s assumption (|4.21|) : 

ll/^+ill —II, ,1+1 11 <- r> lb ,™+ 1 ll 1_<5 ll/i"+ 1 ll (5 

h 

< C s \\u n+1 \\\~ & \\VLj n+l \\ & 2 < C S C^- S C( S . (4.83) 

For simplicity, by taking 8 — |, we see that (|4.21[) is also valid at time step t n+1 if 
we set 



Ci=C S y/C 6 C w . (4.84) 

Note that Ce and Cis are independent of Ci in the derivation. The constant C\ is 
only used in the time step constraint (|4.40p . Therefore, an induction can be applied 
so that the a-priori H 5 assumption f)4 . 2 1 [) is valid at any time step under a global 
time step constraint 

Ati lclkcT a - <4 - 85) 

Again, note that both Cq and Cis contain an exponential decay in time and therefore 
are bounded by a given constant in time. 
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In other words, under (|4.85|) . a global in time constant constraint for the time 
step, the proposed semi-implicit scheme (|4.10p - (|4.12l) is unconditionally stable (in 
terms of spatial grid size and final time). In addition, an asymptotic decay for the L 2 
and H 1 norm for the vorticity (equivalent to H 1 and H 2 norms for the velocity) can 
be derived. Lemma 14.21 is proven. 

Appendix A. A Wente type estimate. The goal here is to present a Wente 
type estimate that is applicable to our doubly periodic setting. Original estimate 
of the Jacobian term (essentially H~ 1 norm) goes back to [45]. Here we need an 
estimate on the L 2 norm of the Jacobian. The case with homogeneous Dirichlet 
boundary condition can be found in [23j [24] . 

Proposition 1. There exists an absolute constant C w > 1 such that 

||VVV0|| H -x <C w \\il>\\ Hl \\<f>\\ H i V^H; er (0),^Ii; er (l]) (A.l) 
||V-h/;-V0|| 2 < C w \m m U\\ m V^H^ r (0),^H; er (Q) (A.2) 

||V x ^-V0|| 2 <C w ||V|| H i||</.|| ff2 G Hp er (0), <f> G Hp er (0). (A.3) 



Proof. Let fi = (0,2tt) 2 as before and ft := (-27r,4vr) 2 . Let p G Cg° (R 2 ) be such 
that p = 1 in O, p = in R 2 — f2 and p(x) € [0, 1] for all x € R 2 . Here ip and <j> are 
27r-periodic functions on R 2 . The proof of (| A.2|) is based on Lemma 1 in [23], which 
states that, in our notation, for pip G Hq(CI) and pep € Hq((1), one has 

||V x (pV) • V(W>)|| L2( ^ < C K (n) IMI^MI^ • (A-4) 

Noting that 

l|v(^)|| ft = ||v(^)|| n + ||vo^)|| fl _ n 

<||V^|| n + ||pVV||^_ n + ||VVp|| 5 _ n 

< WWa + \m\n-u + IMIn-Jl v PlU ( n_n) 

< HWHn + 8 llWHn + 8 ^||V^| fi _ n ||Vp|| iDO(fi _ n) , 



(A.5) 



and a similar computation for ||j0^|| ff2 , the right-hand side of (|A.4j) is majorised as 

WP^W^^Wp^min) ^ (9 + 8col|Vp|| io . (R2) ) 2 C^(^) 2 ||^|| ff2(o) ||0|| H1( ^. (A.6) 

Since the left-hand side of (|A.4p majorises ||V ip- V^|| i2 mp (IA.1I) follows. The proof 
of (| A.3|) is completely analogous, using the estimate [331 Lemma 1], 

||V x (pV) • V(^)|| £a(fl) < C K 11/^11^1(0)11^11^(6) • (A-7) 

for G Hq{&) and G i?£(f2). 

For (jA.ip . we take w G ^o(^) and compute 

HV^-V^l^x^j = sup (V x V-V0,«;) i2(n) 

IMIffl(f2) = l 

< sup ||V0|| £a(fi )||V^| za( n)||VHIi*(n) ( A -8) 

IMIffl(t2) = l 
= l|V0|| L 2 ( ^ ) ||VV'|| i 2 ( ^ ) 
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where the inequality follows from (3.8) in [35]. Arguing as above, (IA.1I) follows. □ 

Appendix B. A convergence result on long time behaviors. Here we 
present a modified version of the abstract result presented in [33] , so that it is appli- 
cable to the current situation, where the phase space is only a subset of a Hilbert (or 
reflexive Banach) space. 

Proposition 2. Let {S(t)} t >o be a continuous semi-group on a complete metric 
space X which is a subset of a separable Hilbert space H with the inherited distance 
(norm) \\-\\. Suppose that the semi-group generates a continuous dissipative dynamical 
system (in the sense of possessing a compact global attractor A) on X . Let {Sk}o<k<k 
be a family of continuous maps on X which generates a family of discrete dissipative 
dynamical system (with global attractor Ak) on X . We further assume that the fol- 
lowing two conditions are satisfied. 

HI: [Uniform boundedness] There exists a k\ G (0,fco] such that {Sk\o<k<k t is 
uniformly bounded in the sense that 

K= |J A k (B.l) 

0<fe<fei 

is bounded in X. 

H2: [Finite time uniform convergence] Sk uniformly converges to S on any finite 
time interval (modulo any initial layer) and uniformly for initial data from 
the global attractor of the scheme in the sense that there exists to > such 
that for any T* > to > 

lim sup \\S%u - S(nk)u\\ = 0. (B.2) 

ueA k ,nke[t ,T*] 

Then the global attractors converge in the sense of Hausdorff semi- distance, i.e. 

lim dist H (A k ,A) = 0. (B.3) 

k— >0 

Moreover, if the following three more stringent conditions are satisfied: 

H3: [Uniform dissipativity] There exists a k\ £ (0,fco) such that {Sk}o<k<k! is 
uniformly dissipative in the sense that 

K= |J A k (B.4) 

0<fc<fci 

is pre-compact in X . 
H4-' [Uniform convergence on the unit time interval] Sk uniformly converges to S 
on the unit time interval (modulo an initial layer) and uniformly for initial 
data from the global attractor of Sk in the sense that for any to £ (0, 1) 

lim sup HSjfcU- S(nk)u\\ = 0. (B.5) 

H5: [Uniform continuity of the continuous system] {S(t)} t >o is uniformly contin- 
uous on K on the unit time interval in the sense that for any T* € [0, 1] 



lim sup ||S(*)u- 5(T*)u|| = 0, (B.6) 
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then the invariant measures of the discrete dynamical system {<Sfc}o<fc<fc converge 
to invariant measures of the continuous dynamical system S. More precisely, let 
fXk € TM.k where IA4k denotes the set of all invariant measures of Sk- There must 
exist a subsequence, still denoted {^h}, and ^ g lAi [an invariant measure of S(t)), 
such that fik weakly converges to fi, i.e., 

/ifc — s- pt, as k — > 0. (B-7) 

Proof. The proof is exactly the same as those in [44j [43] . We leave the detail to 
the interested reader. □ 
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